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List  of  symbols 


Symbols  listed  in  Part  I  of  this  report  are  not  repeated  here. 


b.., . half  base-length  of  wing 

1 . wing  chord 


l0...*...«  distance  between  base  and  x-axis. 
s  ........  exponent  in  representation  of  wing 


cross-section. 

P  .  angle  at  base  of  triangle,  p  =  -  a. 

p(C> . complex  function  replacing  <p(C)  of  Part  I 


General  Remarks 


1.)  Introduction,  In  Part  I  of  this  report1 ^  the  basic 
equations  for  shallow  conical  shells  have  been  derived 
assuming  the  temperature  distribution  in  the  shell  to  be 
symmetric  with  respect  to  the  vertical  plane  of  symmetry 
of  the  delta  wing.  Boundary  conditions  have  been  formulated 
and  a  perturbation  method  has  been  developed  consisting  of 
two  steps  termed  first  and  second  approximation,  respectively. 
Temperature  distribution  has  been  resolved  into  two  parts  5 
one  symmetric  with  respect  to  the  horizontal  plane  of 
symmetry  of  the  delta  wing,  and  the  other  skew-symmetric 
with  respect  to  that  plane.  The  general  solution  and  a 
numerical  example  have  been  given  for  the  symmetric  part 
using  first  approximation  only. 

In  the  present  Part  II  of  the  report  the  study  of 
thermal  stresses  in  delta  wings  is  continued.  Equations  of 
first  approximation  for  the  skew-symmetric  temperature  part 
are  derived  and  solved.  A  numerical  example  is  given.  Then 
the  equations  of  second  (and  final)  approximation  are 
obtained  and  solved  for  both  the  symmetric  and  the  skew- 
symmetric  part.  The  numerical  example  is  continued  into 
the  second  approximation. 

Frequent  reference  is  made  in  the  following  to  Part  I. 
Equations  or  figures  contained  in  Part  I  are  quoted  by 
adding  a  Roman  I  to  the  equation  number  or  figure  numbers 
Eq. (1-12-1)  or  Fig. 1-3. 


1) 


Parkus  (see  list  of  references) 


Skew-Symmetric  Temperature  Distribution 


2.)  General  Equations.  The  two  basic  equations  (1-7-2) 
remain,  of  course,  valid  in  the  present  case  of  a 
temperature  distribution  skew-symmetric  with  respect  to 
the  x,y-plane,  Fig.1.  Boundary  conditions  (1-8-3),  however, 
have  to  be  changed.  Since  upper  and  lower  shell  now 
experience  the  same  deflection  w  but  opposite  tangential 
displacements  u  and  v,  Fig. 2,  the  conditions 


'i 


> 


u  =  0 


v  =  0 


(2-1) 


have  to  be  satisfied  along  base  (upper  sign)  and  legs 
(lower  sign)  of  the  triangle,  cf. Fig.1.  Boundary  conditions 
(1-8-4)  along  the  y-axis  (line  of  symmetry)  remain 
unchanged : 


dw 

3x  = 


0 


\ 

i 

| 

u  =  0 


nxy  =  0 


■s 

j 
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(2-2) 


Equations  (2-1)  have  now  to  be  expressed  in  terms  of 
the  deflection  w  and  the  stress  function  F.  Eqs. (1-7-5) 
and  (1-7-6)  render 


d 


2 


Furthermore 


D 


q  =  -  I  £-(v2w  +  a'  ) 

9  r  1 


Base  and  leg  of  the  triangle  are  given  by  9  -  0  and 
9  =  j?  =  ^  -  a,  respectively.  Hence, the  first  two  of 
conditions  (2-1)  take  on  the  form 


(2-4) 


V2  W  —  (  1  —  V )  ——IT  +  AT.  =  0 

dr^  1 


i-fv2.  +  nr  )  ♦  (1-V)r  &)  -  1  ^  ^  =  0 

69  Ar  Ac p  dtp  dr 


>  (2-5a) 


along  cp  =  0  and  9  =  j!  . 

From  eqs.( 1-5-3)  and  (1-3-7)  one  finds  using  polar 
coordinates  and  denoting  the  radial  displacement  by  u  and 


d2F 

d7 


the  change  of 

du 

1  \r£l* 

dr 

®LVF 

—  + 

^  +  £  , 

r 

59  r 

+ 

r2  de  = 

d<f 

dr 

Along 

0 

11 

0- 

EH 


dr‘ 


>  (2-6) 


^  V"  '  * 


all  values  of  r.  Hence,  from  the  first  equation 

2 

V2F  -  (1+  v)  £-|  +  ( 1  -v)B  %o  =  0 
dr 


(2-5b) 
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Melan-Parkus,  p.59 
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Differentiating  the  first  equation  with  respect  to  9  and 
the  third  with  respect  to  r,  and  subtracting,  there  follows 

—  Sv2P+  (1+v)r  £-*-(|)  +  (l-v)Da*  =0  (2-5c) 

d<p  l  dr 

along  9=0  and  9  =  j5. 

The  four  equations  (2-5a)  to  (2-5c)  replace  the  boundary- 
conditions  (2-1). 


3.)  Perturbation  procedure.  Expanding  w  and  P  in  series 
in  powers  of  c,  Eqs.( 1-9-2),  and  retaining  only  the  first 
two  terms  one  obtains  eqs. (1-9-3)  and  (1-9-4)* 


2  2  p  2 
V  =  f  -  V^1 

V2V2P1  =  -  (  1-v)DV2  X 


(3-1) 


2  2 

IV  7  w. 


V2V2Pn 


1  1,2  h 


dr 


Eh 

r  dr2 


(3-2) 


Performing  the  same  expansion  in  the  boundary  conditions 
(2-5a)  to  (2-5c)  and  using  ( 1—9—1 )  one  finds 
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2  ft2pl 

y^P1  -  (1+v)  — gi  +  (l-v)D  X0  »  0 
dr 


(3-0 


(3-5) 


(3-6) 


The  four  pairs  of  boundary  conditions  (3-3)  to  (3-6)  are 
valid  along  base  and  legs  of  the  triangle,  i.e.  for  <p  =  0 
and  <p  =  (3.  The  remaining  four  conditions  (2-2),  valid  along 
the  line  of  symmetry  x  =  0,  take  on  the  form  (1-9-7)  and 
(1-9-8)  and,  as  in  the  case  of  the  symmetric  temperature 
distribution,  can  be  taken  care  of  automatically  by  extending 
the  validity  of  the  preceding  equations  to  the  entire 
triangular  region  of  the  wing.  For  the  second  approximation, 
however,  an  additional  external  load 
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4.)  First  approximation.  The  solutions  of  eqs. (3-1)  are 
split  up  into  particular  solutions  F^  of  the  complete 
equations  and  solutions  w^,  F^  of  the  homogeneous  equations 
If,  as  in  Part  I  of  this  report,  attention  is  restricted  to 
temperature  changes  only  then  p  =  0  in  eq.(3-1)  and  the 
particular  solutions  may  be  obtained  from  eqs. (1-10-2) : 

\ 


V2F^  =  -  (l-v)D  * 


(4-1) 


No  boundary  conditions  need  be  taken  into  consideration. 

Eqs. (4-1)  can  be  solved  by  the  same  method  as  used  in 
sec. 10  of  Part  I.  Another  possibility  consists  in  employing 
the  methodof  finite  differences. 

The  particular  solutions  thus  obtained  will,  in  general, 
not  satisfy  the  boundary  conditions  (3-3)  and  (3-4). 
Substituting  =  wtj  +  into  (3-3)  and  making  use  of 
eqs. (4-1)  one  obtains 


The  same  functional  representation  is  assumed  for 
as  giver,  by  eq.(l-11-3) 


(4-2) 


2w^  =  zp(z)  +  zp(z)  +  ^(z)  +%(*) 


(4-3) 
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where  p(z)  and  ^  (z)  are  analytic  functions1 ^  of  the  complex 
variable 

z  =  x  +  iy  =  re1'’5  -  a,  a  =  1  tan  a  +  ilQ  (4-4) 

A  bar  denotes  the  conjugate  complex  quantity.  Using 

2-  =  e1?  i-  +  e"i{P  1-  ,  11  =  i(ei<P  1.  _  e*i<P  1.) 

<5r  dz  dz  r  d<p  dz  dz 

and  writing 

%  »(z)  =y  (z)  (4-5) 

one  finds 

v2w^  =  2  p'(z)  +  p'(z)  (4-6) 

1  L.  V2w^  =  21  p"  (z)ei(P  -  p"Yz)e~i<? 

±  .  dw^  r  It 

— <r  dV^  =  ?  ^  (z)  +  e2l<P  -  >  (4-7) 

-  |  zp"  (z)  +  y  '(z)  e"2iiP 

^2wh  _  _ ^ 

jpj-  s  £  zp"(z)+Y’(z)  e2i<t>  +  £  zp"  (z)  +  y'(z)  e*fc2i(,p  + 

_  (4-8) 

+  p'(z)  +  p* (z) 

1 ) 

The  letter  p  is  now  used  instead  of  <p  to  avoid 
confusion  with  the  polar  angle. 
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In  order  to  simplify  the  boundary  conditions  the  second 
of  eqs.(4-2)  is  multiplied  by  dr  =  e^^dz  =  e^dz,  and 
integrated: 


2i 


ji'(z)  -  p.'  ( z)  +  j  (1-v)<j  zp”(z)  +  Y'(z) 


32i?  - 


zp”  (z)  +  Y'(z) 


e'21*r  -  (i-v,fe$ +  o'(^. 


is  real*  The  first  of  eqs*(4-2)  assumes  the  form 


(1+v)  jp'(z)  +  p'(z)j  -  zp”  (z)  +  y'(z) 


32i!p  + 


+  zp"  (z)  +  Y  '  (z) 


,-219  l  = 


d2wP 
=  (1-v)  —4 
dr 


Multiplying  this  equation  by  i  and  adding  it  to  equation 
above  renders 


C5+v)p'(z)  -  (l-v)p'(z)  -  (1-v)  zp"(z)  +y'(z) 

=  (1-v)  —J-  +  i  — (rr1)  +  i  0.(9) 

dr^  drd9  r  1 


,-2if 


This  equation  is  now  once  more  integrated  with  respect  to  r. 
One  has  with  eq.(4-4) 

Jzp"  (z) e""2i<?dr  =-J{%  +  a  -  ae“2i^)ji!'  (z)dr  = 

*  e^J  zp”(z)dz  +  (aei(?  -  aei^)J'  p”  (z)dz 

=  ei<p  zp'(z)  -  p(z)~|  +  (aei<p  -  ae^^Jp’ (zjae^^zp' (z)  -  e^pf*) 


and  hence 


X)i(z)  -  zji'(z)  -y(z) 
where 


=  e1? 


dw^ 

dr 


+  i 


A  W1 

— (— ) 

dip r  _ 


i  C.j((p)z  +  C2(q») 


(4-9) 


The  triangular  region  is  simply  connected  and  and  02  may 
he  put  equal  to  zero.  Therefore  the  boundary  condition  reads 

finally  __  _  j„P  d„P 

-y(z)  =  (-1  +  i  1  — i)ei(P  (4-10) 


a,p(z)  -  zp'(z)  -  y(! 


dr 


dip 


along  9=0  and  <p  = 

For  the  stress  function  both  differential  equation 
and  boundary  conditions  are  the  same  as  for  the  deflection  w^ . 
Hence,  upon  putting  as  in  eq.(l-1 1 —1 9 ) , 


2F *  =  z$(z)  +  zj(z)  +  X(z)  +  X(z) 


(4-11) 


i/(z)  *  X* (z) 


one  has  at  once  from  eqs.(3-4) 


<JF?  <  dF? 


*,|(z)  -  z{*(z)  -1J/(Z)  =  (*~-  +  i  |  (4-12) 


where 


-  10  - 


5.  Example.  Consider  the  particular  case  where  the 

temperature  increases  from  its  initial  uniform  value  zero 

to  a  constant  value  TQ  in  the  upper  half  of  the  wing  and, 

at  the  same  time,  decreases  to  the  value  -T  in  the  lower 

o 

half.  Then 

r0  -  (1+v)aT0,  T1  =  0 

and,  from  eqs.(4-1),  taking  symmetry  with  respect  to 
x  =  0  into  account, 

=  0,  =  -  §(x2+y2)  =  -  £  zz  (5-1) 


where 


(5-2) 


Substitution  into  eq.(4-10)  renders  w^  =  0  and  hence 
w1  =  0,  while  from  eq.(4-12)  one  has 


t*$(z)  -  z$'(z)  -  Yz)  =  “  Cz 


(5-3) 


The  solution  of  this  equation  can  easily  be  given  in  closed 
forms 

$(z)  =  “flPT  z»  *  0  (5-4) 

Eq.(5-3)  is  then  satisfied  not  only  on  the  boundary  of  the 
triangular  region  but  everywhere  within  this  region.  Prom 
eq.(4-11)  one  finds 


*h  „  C  - 
*1  ■  zz 


(5-5) 


11  *- 


and  hence,  using  eq.(5-1), 

F.j  =  “  ^  zz  ( 5-6) 


The  corresponding  state  of  stress  is  uniform  compression 
in  the  upper  half  and  uniform  tension  in  the  lower  half  of 
the  wing.  In  particular,  one  finds  in  the  upper  shell 


d2F„ 


nx1  = 


ay 


tf+1  n  _  „ 
'«=T  C  '  nyl  * 


nxyi  =0 


(5-7) 


In  the  case  of  a  more  complicated  temperature  distribution 
no  closed  solution  can,  in  general,  be  found  and  a  numerical 
procedure  has  to.  be  used.  In  order  to  demonstrate  this  the 
problem  stated  above  will  now  be  solved  numerically. 

Eq.(5-3)  is  transformed  from  the  triangular  region  into 
the  circular  domain  by  means  of  the  mapping  function 
z  =  co(£),  cf.eq. (1-11-6), 


x}(6T)  -  -'UKff)  =  CH(60  =  -  Coa(er)  (5-8) 

u)'(ff)  1 

valid  on  the  boundary  £  =  Ef  of  the  unit  circle.  Using  the 
same  expansions  as  in  eqs . (1-1 1-8a)  and  (1-11-22)  to 
(1-11-24)  and  utilizing  the  Cauchy  formula  (1-11-12)  the 
following  equation  is  obtained 

°o  00  rtf  1  _  00 

xYlk)^k  ~  YZ  Z _ m  \  °mfk-1  ^  ”  ®o  *  °^Z]Hkck 

1  k=0  m=1  0 

Comparison  of  coefficients  of  equal  powers  of  £  renders 
the  following  set  of  linear  equations  for  the  expansion 
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coefficients  Ak  of  funotion  |(£) 

OqI,  +  20^X2  +  ...  +  (n+1)on  +  Eq  *  -  Hq0  (5-9) 


( 


end  n 

*Ak  -  J^m  cmfk_1  ^  »  HkC  (k  -  1,2,..n)  (5-10) 

m=1 


A  similar  set  of  equations  for  the  coefficients  ^  of  the 
function  ^  (£)  is  obtained  by  replacing  eq.(5-8)  by  its 
conjugate  complex  and  performing  the  same  operations  as 
before.  One  gets 


n 


c!U  + 


(k=1,2...2n-1)  (5-11) 


Coefficients  on  are  all  real  quantities  and  have  already 
been  tabulated  in  Part  I,  Table  12,  for  a  set  of  10  triangles. 
Coefficients  Hk  of  the  expansion  (1-11-22) 

+n 

H(fl’)  HkSk  (5-12) 

-n 

have  to  be  calculated  by  means  of  Fourier  expansion  of  H(C) 
along  the  unit  circle,  cf. sec. 1-16.  In  the  present  example 
we  have 

H(S)  =  -u>(ff)  -  -  ReiY  -  -  R(cosY  +isinY) 


Real  and  imaginary  part  of  H( CT)  as  functions  of  the  polar 
angle  &  on  the  circumference  of  the  unit  circle  can  be 
taken  from  tables  1-1  to  1-10,  of.  Fig. 1-10.  As  in  the 
example  in  Part  I  a  triangle  with  angle  2a  *  100°  is  assumed. 


.1 


Pigs  .3  and  4  show  the  two  functions  -RejH(ff)^  and 
-  Im  £(*)}.  A  Pourier  analysis  renders  the  corresponding  Pourier 
coefficients  which  in  turn  determine  the  coefficients  H^,  by 
means  of  the  equations 


Ho  -  *V  \  H-n  *  ?<Vbn)  <5-15> 

The  Pourier  coefficients  are  given  on  Pigs. 3  and  4  for  n  *  6 

while  the  coefficients  are  contained  in  Table  1 .  The  approximation 

curves  are  also  represented  (broken  lines). 


k=0 

1 

2 

5 

4 

5 

e 

p 

Hk 

-0,26i 

-10,30i 

1,80i 

0,43i 

i 

vo 

H* 

0,68i 

0,191 

0,047671 

H-k 

0 

-0,08i 

0,07i 

0,02i 

-0,081 

0,07i 

Table  1 .  Power  series  coefficients  of  H 
Each  value  in  the  tables  1  and  2  is  to  be  multiplied  by  the 
corresponding  factor  ji. 

The  exact  value  of  all  quantities  is  zero.  Table  1 
reflects  the  error  involved  in  the  approximation  of  the  mapping 
function. 

Eqs*(5-9)»  (5-10)  and  (5-11)  can  now  be  set  up  and  solved. 
The  results  are  given  in  Table  2.  One  notices  that  the  are 
not  precisely  proportional  to  the  and  the  are  not  all 
zero  as  would  be  the  oase  with  the  exact  solution.  The  error 
is  due  to  the  relatively  small  dumber  n  *  6  of  terms  in  the 
series  expansion. 
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k=0 

1 

2 

3 

4 

5 

6 

m 

□ 

-8,874i 

1,417i 

0 , 308i 

-1,334i 

mm 

*k 

-0,770i 

-0,677i 

2,086i 

-0,1 17i 

-1,369i 

4,818i 

Pllll 

0,0476701 

■ 

| 

8 

9 

10 

11 

Bk 

-0,579i 

0,356i 

-0,076i 

-0,1 07i 

-0,012i 

Table  2.  Power  series  coefficients  of  $  and 


( 


Second  Approximation 

6.)  Symmetric  temperature  distribution.  As  can  be  seen  from 
a  comparison  of  eqs. (1-9-4)  and  the  corresponding  boundary 
conditions  (1-9-6)  with  eqs. (1-9-3)  and  boundary  conditions 
(1-9-5)  the  steps  to  be  performed  in  the  second  approximation 
are,  in  principle,  the  same  as  in  the  first  approximation.  They 
deviate,  however,  in  two  respects.  First,  the  equations  for 
the  particular  solutions  and  5^  can  no  longer  be  reduced 
to  the  simple  form  (1-10-2).  Second,  a  line  load  has  to  be  placed 
along  the  axis  of  y  in  order  to  have  the  boundary  conditions 
satisfied  there. 

The  determination  of  the  particular  solutions 


o_2  n  1  dZpi  d2,1 

Nv  Vw£  =  -  1  (— -1  +  — wi) 

*  T  dr*  dr * 

2  2  n  Bh  ^W1  d2w£ 


(6-1) 


may  be  performed  either  in  the  triangular  s-plane  or  in  the 


circular  £-plane,  where  z  =  w(£) .  Since  and  are  uaually 
known  in  the  z-plane  a  series  expansion  similar  to  (1-10-5) 
and  (1-10-7)  may  be  used  for  the  corresponding  parts  and 


5*2  -|  of  w|  and  igt 


W?1  =  n  008 

2 1  'Bn 


nittp 


p?,  =  p.  ="*k 


‘21  mn 


and,  in  case  where  m  +  k  =  —  , 
1  a 


cos 


nit<p 


(6-2) 


w|i  =  ^r111*11  In  ^  cos(m+k)<p 
=  5mrmfk  In  ^  cos(nH-k)9 


(6-3) 


In  contradistinction,  and  pj  are  known  in  the  £-plane . 

The  corresponding  parts  w|2  and  F^  w2  may  "therefore 

more  conveniently  be  represented  by  series  expansions  in  terms 
of  9  and  k  .  Using  the  relations 

A  =  0(9, 


V2  =4  £  -  — 4 -  ± 


A  = 


tzbi  ut'io^Tcj  dZdl  w*  (Om’TST 


=  %+iL+J_i_,  e(o,^)  - - — 

To*  ?  dp  0^ 


i o'U)&Tc ) 


(6-4) 


and  taking  eq,.(4-8)  into  account  eq,sj(6-1)  may  be  written 


Na(©Aw|2)  <=  Rfrj 


«T(C>  1 


e21^  +  2f  (C) 


(6-5) 
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a(0^2|2)  *  H  Re^ 


^51  f'U)  +  *?*(£)  e2i*  +  2fU) 
CO'  (C)  r 


>  (6-6) 


where,  in  accordance  with  the  notation  of  eqs. (1-1 1-28), 


F(«)  =  ,  p*(o  -  y*(o  -  ~ 

w'(0  I  W’(S)  r  u)'(C)  «*>' 


f'(C) 


(« 


C 


In  order  to  evaluate  the  right-hahd  side  of  eqs.  (6^*5)  and 
(6-6)  the  function  l/w'(£)  has  to  be  expanded  into  a  series 


O)'(S) 


=L»„«n-E«n?n«ln*  (6-7) 


[The  coefficients  an  in  this  expansion  may  be  found  by  putting 
p  =  1  and  using  eq.( 1-13-2)  which  gives  the  values  of  c O' 
along  the  circumference  of  the  unit  circle.  From  eq.(6-4)  one 
has  then 

eCj.JO  “ILwn+mei(,1'B)'5'=IEPr,  ?r°°»  (6-8) 

r  s 

The  solution  of  eqs. (6-5)  and  (6-6)  proceeds  now  in  the 
following  manner.  The  numerical  values  of  the  right-hand  side* 
of  these  equations  are  calculated  at  a  sufficient  number  of 
points  in  the  C-plane  using  a  po?.ar  coordinate  network  ft/fr  • 
From  these  values  a  series  expansion  is  obtained  with  the  aid 
of  some  suitable  numerical  method: 


J  f  f 1 


2h  R 
r9  Re<i 


.*>'(€) 

■«(o 

cTic) 


(o  +y*(s) 


e2i*  +  2$*(£)  f  -  Hwm~2<s°8 


}**'(£)  +Y*(ol  e2i<p  +  2j»*(0  L  »  ^Hamn?m”2°08 


X6-9) 
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( 


The  expansion  will  contain  cosine  terns  only  sinoe  all  functions 
are  symmetric  with  respect  to  the  y-axis.  A  similar  expansion  ■» 
with  unknown  coefficients  -  is  then  assumed  for  the  following 
quantities 


0ow|2  83  X^'bmnPmcos  n^+  ZZ  ^n?11  ln  P  008 

m  n  n  [  (6-10) 


6n5f2  =  2Z^BnmPmcos  n^+  H  Bn9*  ln  P  008  J 
m  n  n  J 


where  the  prime  at  the  summation  sign  indioates  that  terms 
with  m  =  n  have  to  be  omitted.  Prom  the  definition  of  the 
operators,  eq. (6-4),  there  follows 


a(0aw|2)  =»XI  YL  (m2-n2)bmn5>m“2  cos  n£  +  2  ZZ  n 


bn5>n“2cos  n/ 


m 


n 


n 


< 


a(0af|2)  a  y  y~ *  ( m2 -n2 ) Bmn9m~2  cos  n£  +  2  ^  n  B^^oos  n/fr 
m  n  n 


Comparing  coefficients  with  eq.(6-9)  one  finds 


mn 


"San 

~~2  ? 

m  -n 


>  (m  ^  n) 


(6-11) 


B, 


amn 


mn  J2r~Z 
m  — n 


B_ 


ann 

2n~ 


In  order  to  obtain  w22  and  1^2  the  foregoing  steps  would 
have  to.  be  repeated.  This,  however,  leads  to  very  time-oonsuming 
calculations.  We  proceed  therefore  in  an  approximate  fashion 
by  applying  the  same  method  used  to  obtain  eq.  (6-9),  i.e*  we 
evaluate  eqs.(6-10)  and  (6-8)  at  all  points  of  the  ooordinate 


18  - 


network  and  expands 

a”22  "51  5jw”"2  008 

iP22  "L^mnP”'2  808  8 A 

Ihen,  upon  putting 

w22  =HL‘wmn?m  008  +  XT^?"111  P  cos 

*22  =F  ZlpmnPm  008  n^+  H  Pn?nin  ?  003  nJ" 


(6-12) 


r (6-13) 


and  substituting  into  eqc.(6-12)  one  finds,  corresponding  to 
eqr .( 6-1 1 ) , 


mn 

wmn  “  ~T~7 
m  -n 


w  = 
n  2n 


> (m  /  n) 


Cr 

m^-n" 


A 


F  = 
n 


nn 

2n“ 


(6-14) 


To  w|  a  second  solution  pertaining  to  the  line  load 
(3-7)  has  to  be  added.  The  solution  for  an  infinite  plate  strip 
simply  supported  along  its  edges  and  carrying  a  line  load  is 
available  in  the  literature^  and  may  be  used  here.  Introducing 
a  coordinate  system  x,y  as  shown  in  Fig.1  one  has 

nuixt  _ * 


q  1'  *n, 

w2  ■  - T  >  -t( 

C  A  W«r  *  L n  ' 


1  + 


4NiT  n' 


niuxt 

UUL).—r-.  tn^-  (6-15) 


T7 


Girkmann,  p.170. 
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where  y*  «  y  +  1_.  The  coefficients  a„  are  those  of  the 
o  n 

Fourier  expansion  of  the  line  load, 


q(y*)  =  -  21 


df 

dx 


n. 


xl 


-r 


8n  sin 


nwy 


(6-16) 


x=0 


For  a  supersonic  wing  profile  the  following  approximation  may 
be  used 


If  =  (5s-  9S)y*  =  (£s-  <p8)r  sin  9  (6-17) 


where  s  is  a  suitable  chosen  positive  number.  Then 


df  a  i  o  6f  „  a  . 

1  —  =  s<p  1  sin^<p,  1  r;  =  ^  -  <p  '(tp  +  f  sin2<p)  (6-18) 
dx  oy  c 


The  ratio  of  the  leading-edge  slope  to  the  trailing-edge  slope, 
taken  along  sections  x  =  konst,  is  given  by  s  sin  2jJ/2|I. 

The  normal  force  n^  at  x  =  0  is  the  sum  of  the  particular 
part  and  of  the  homogeneous  part  n^.  The  first  is  in 
terms  of  F^,  eq..  (1-10-5) ,  given  by 


9  =  0 


while  for  the  second  one  finds  from  eqs. (1-1 1-30) 


(6-19) 


(6-20) 


The  coefficients  an  in  the  expansion  (6-16)  may  now  be  calculated 
by  some  suitable  numerical  method. 

The  homogeneous  solutions  w|  and  Pg  of  the  seoond  approximation 


-  20  *- 


are  introduced  in  the  sane  form  as  in  the  first 
approximations,  eqs.(I-11-3)  and  (1-11-19) 


2w|  =  z  ji( z )  +  zjKz)  +^(z)  +  (z) 


( 6-2 1 ) 


2F^  =  z$(z)  +  z^(z)  +  X(z)  +  X(z)j 


The  corresponding  boundary  conditions,  eqa(l-9-6),  are 
identical  with  those  of  the  first  approximation,  eqa(  1-9-5). 
Hence,  eq.. ( 1-1 1-4)  remains  unchanged: 

_  _  ^  ^  9. 

u(z)  +  zu‘(z)  +  -y(z)  =  -  ( — -  +  i  — -)  (6-22) 

r  r  d*  by 

or,  after  transformation  into  the  £-plane, 

^(trs  -  -  i) w9  P4-11 

u(ef)  +  n'(er)  +  *(*>  -  -  (—  +  i  — )  =  &«0  (6-23) 

co’(er)  dx  dy 

Similarly,  for 

{(ff)  +  JiLiSl  $•(*)  +^(5)  =  -  (^+  i  ^k)P-  H(fO  (6-24) 

«•(*)  1  bx  fly 

We  note  that 


a 

3x 


+  i 


fc  - *lT<* 


+  i  a  \ 

3r  +  r 


(6-25) 


7.  Skew-symmetric  temperature  distribution.  2qs.(3-2)  are 
valid  in  this  case.  They  are  identical  with  those  of  the 
symmetric  distribution.  Moreover,  the  same  line  load  as 
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( 


( 


in  the  aymmetrio  ease  given  by  eq.(6-l6)  has  to  be  placed 
on  the  axis  of  y.  Boundary  conditions,  however,  are  now  In 
the  form  of  eqs.(3-5)  and  (3-6). 

The  particular  integrals  w|  +  w|  and  are 

determined  in  exactly  the  same  way  as  in  the  symmetric  case, 


sec. 6. 


Por  the  homogeneous  solution  w£  one  has  from  boundary 
donditions  (3-5)  on  qp  =*  0  and  <p  =  j5  ,  with  w|  +  w|  *  w*  , 


V^Wg  -  (  1-v)  — yp  =  (1-v) 
*  dr^ 


d2w* 

0  2  „2  • 

-  w  -  V  W0 

dTd  £ 


d 

dtp  L 


1  .  (1-v) 

1  C  i>T  37 


>  (7-1) 


1(1  0,  l!fl  .  4.  f  V.! 

“  r  d<p  dr2  Lr 


2  *  (1-V) 


n 


teg  is  now  split  up  into  two  parts, 


_  ** 
"2  "21 


+  w* 
"22 


(7-2) 


where 


2wJ1  =  2)1.,  (2)  +  2^(2)  +Xi^z)  +  Xl(z>  (7-3) 

2w22  =  2P2(z)  +  ■fef*)  +  J('2^  +  XV8)  (7-4) 

The  first  part  w^  is  obtained  as  a  biharmonio  funotion 
satisfying  boundary  conditions  (4-2)  or  (4-10),  with  w^ 
replaced  by  w*  .  The  second  part  w^g  is  then  subject  to  the 
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( 


following  boundary  conditions 


*2  h 

-  d-v) 


dr 


-  v2*^ 


d_ 

df 


1  1  df  d^F, 
-  ~>  — 1 


df  dr 


-  v2**) 


dip 


P (7-5) 


We  proceed  now  as  in  see. 4.  The  second  of  eqs.(7-5)  ie  multiplied 
by  dr  and  integrated.  Bearing  in  mind  that  df/rdf  is 
independent  of  r,  of.  eq.(6-17),  and  putting  the  integration 
constant  equal  to  zero,  we  obtain 


Jr. 


2i|p£(z)  -  }i£(z)J  +  |(1_v)^[zp2(z)  +yiM]  e2i(p 

[•Sw  *^>] °'21]-  W  -J^<v2,*,ar 


This  equation  is  multiplied  by  i  and  then  subtracted  from 
the  first  of  eqs.(7-5),  to  yield 


Ov)p^(z)  -  (l-v)ji^(z)  -  (1-v)  zp^(z) 


Another  integration  with  respeot  to  r  renders 
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*P2<b>  ’  ZP’<Z>  -Y2(z)  ■  *  2Jv2w2dr]f=^  (7“6) 

where  \  is  given  by  eq.(4-9). 

The  same  procedure  may  be  applied  to  boundary  condition 
(3-6).  letting 

*2  *  *21  +  *22  (7-7) 


and  putting 

( 

2I^1  «  z^(z)  +  z^(z)  +  X.,(z)  +  X.j  (z) 

(7-8) 

2F^2  »  z|2(z)  +  z|2(z)  +  X2(z)  +  Xg(z) 

(7-9) 

one  has  boundary  condition  (4-12),  with  F^  replaced  by 
f|,  valid  for  while  F^2  ha8  to  satisfy  the  following 

equation 

tt<f2(z)  -  z$£(z)  -^2(z)  =  -  2 

where  3C  is  given  by  eq.(4-l3). 

8,  Example.  The  example  of  sec. 5  shall  now  be  continued  into 
the  second  approximation.  With  w1  *  0  one  finds  at  once  from 
the  second  of  eqs.(3-2)  together  with  boundary  conditions  (3-6) 

F2  =  0  (8-1) 

In  order  to  determine  w2  one  has,  however,  'to  perform  all  the 
steps  outlined  in  sec. 7.  First  the  particular  solution  w| 
must  be  found.  Substituting  eq,(5-6)  into  the  first  of  eqs. 
(6-1)  one  readily  gets 


I' 


72F?dr 


e1? 

THJ 


(7-  10) 
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w2  “  rjfl  grf  +  (8-2) 

where  B(r)  is  some  arbitrary  biharmonic  function  which  will 
be  determined  later. 

Next,  the  solution  w|  due  to  the  line  load  (3-7)  has  to 
be  calculated.  Substituting  nx1  from  eq.(5-7)  into  eq.(6-16) 
and  using  the  first  of  eqs.(6-18)  one  obtains 

q(y*)  =  2C  sqps—^ sin2<p  (8-3) 

( 

A  value  of  s  =  0.08  has  been  chosen  in  the  following.  Pig. 5 
shows  the  corresponding  cross  section  of  the  wing  along 
x  =  const.  Two  other  types  with  s  =  0.07  and  0.1,  respectively, 
are  also  shown.  The  relative  maximum  thickness  of  the  wing 
in  the  three  cases  is  4.0  $,  4*6  #  and  5.6  respectively. 

In  accordance  with  eq.(6-l6)  expression  (8-3)  is  now 
developed  into  a  Fourier  sine  series.  As  has  been  done 
consequently  in  this  example  only  the  first  6  terms  in  the 
series  are  retained.  Table  3  shows  the  corresponding  coefficients 


ak’  wi'fch  v  =  0.3  in  the  expression  (qe +1  )/(* -1)  s  2/(1-v). 


k 

1 

2 

3 

4 

5 

6 

V° 

0.1959 

-0.0868 

0.0556 

-0.0426 

0.0333 

-0.0280 

Table  3. 


For  the  calculation  of  the  first  part  wj^  of  the  homogeneous 
solution  Wg  eq.(4-10)  is  used  in  the  form 

-  -  dwt  dwt  / .  \ 

Ajl.j(z)  -  ZJlj(z)  -  ^(z)  =  35-  +  i  -gj-  »  G^Cff) 


(8-4) 


dw^/dx 

dw|/dy 

dw|/dx 

dw|/dy 

dw*/dx 

dw*/dy 

n 

181 .87 

152.62 

mom 

181.87 

150.07 

m 

95.19 

79.69 

0.18 

-2.29 

95.57 

77.40 

15 

59.49 

49.72 

0.51 

-2.06 

59.80 

47.66 

50 

28.48 

25.82 

0.48 

-1.69 

28.96 

22.15 

45 

4.65 

5.90 

0.58 

-1 .58 

5.25 

2.52 

60 

-15.85 

-15.27 

0.66 

-1.08 

-15.17 

-14.35 

75 

-55.29 

-29.54 

0.71 

-0.78 

-54.58 

-30.32 

90 

-52.94 

-44.41 

0.72 

-0.46 

-52.22 

-44.87 

105 

-68.51 

-57.66 

0.67 

-0.11 

-67.84 

-57.77 

112.5 

-89.92 

-75.28 

0.58 

0.1.2 

-89.54 

-75.16 

120 

0 

0 

0 

0.24 

0 

0.24 

127.5 

-116.44 

0 

0 

1.15 

-116.44 

1.15 

155 

-88.68 

0 

0 

1.45 

-  88.68 

1.43 

150 

-64.91 

0 

0 

1.75 

-  64 . 91 

1.75 

165 

-54.77 

0 

0 

1*95 

-  54.77 

1.92 

160 

0 

0 

0 

1  *99 

0 

1.99 

Table  4. 

All  values  are  to  be  multiplied  by  C/N. 
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with  w*  =  w|  +  w^. 

Table  4  contains  the  values  of  the  derivatives  of  Nw|/C 
from  eq.(8-2),  of  NwJf/C  from  eq.(6-15)  and  of  Nw*/C  as 
functions  of  the  polar  angle  &  in  the  5-plane,  see  Fig.I-*6. 

The  as  yet  undetermined  function  B(r)  in  eq,(8-2)  has  now  been 
chosen  in  such  a  way  as  to  make  the  discontinuity  produced 
by  P  in  ^  =  180°  vanish.  One  finds  B(r)  *  “  jjl  +■}  br2. 

The  corresponding  graph  is  represented  in  full  lines  in  Figs. 6 
and  7. 

Real  and  imaginary  part  of  G^(5)  in  eq.(8-4)  are  now 
developed  into  Fourier  series,  represented  by  the  broken  lines 
in  Figs. 6  and  7.  The  coefficients  Gn  in  the  expansion.  (1-11-9) 
follow  then  from  the  equations 


41’  'lao<  “ij’  "  2< W  <8'5) 


They  are  all  purely  imaginary  and  are  given  in  Table  5. 
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Eqs.(8-4)  are  now  solved  by  introducing  expansions  for  ^  and 
y1  of  the  form  (1-11-10)  and  (1-11-11).  The  results  are 
presented  in  Table  6. 


n 

0 

1 

2 

5 

4 

5 

6 

a*1* 

n 

0 

-  9.97 

1.24 

0.27 

1.03 

0.32 

0.39 

b<’> 

n 

-1.86 

-17.96 

59.60 

11.30 

8.40 

16.88 

5.73 

( 


r. 

7 

8 

9 

10 

11 

n 

-1 .98 

-  0.82 

-0.36 

-0.28 

-0.07 

Table  6. 

All  values  are  to  be  multiplied  by  i  ^ 

For  the  calculation  of  the  second  part  of  the  homogeneous 
solution  the  foregoing  procedure  has  to  be  applied  to  eq.(7-6). 
TJsing  eq.(6-15)  one  finds 

_  _  nitjxl  g 

^w2  *  “  dfet  n"  e  8in  tF” 


This  expression  is  zero  along  the  base  y  *  -  1Q  of  the 
triangle  while  along  the  legs  one  has 


2 

1  cos 


i  \  r.* 

n 


2Nit 


(oos 


_  tan  a  sin  £^JL_) 


-  e 


-nit tan  a 


£q.(8-2)  renders 


„2mp  C  afcj-1/  2  .  v 
^w2  "  IT  >r=T(r  "  7  b) 

whence 

r 

T  2  p,  C  *!+t  ,r  2  vv 

J  v  w2dr  =  IT  r(?  “  3  b) 

o 

In  addition,  from  eq..(6-17), 

|  =  (ps-<ps)cos  cp  -  s?8""1  sin  9 

Using  eq.(5-6)  for  F^  the  right-hand  side  of  eq.(7-6),  denoted 
hy  may  now  be  calculated.  Pigs. 8  and  9  show  the  graphs 

of  its  real  and  imaginary  parts  together  with  the  corresponding 
Fourier  expansions  (broken  lines).  Prom  these  coefficients 

are  found  with  the  aid  of  eqs,(8-5)  and  are  represented  in 
Table  7. 


k 

-6 

-5 

-4 

-5 

-2 

-1 

42> 

-9.0 

62.1 

.  .  _ 

-70.2 

175.0 

-200.1 

556.0 

k 

B 

1 

2 

5 

4 

5 

6 

SI 

457.6 

-454.5 

-25.7 

-159.2 

65.4 

-125.1 

17.4 

Table  7 


All  values  are  to  be  multiplied  by  i.C/S, 
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The  last  step  consists  in  the  determination  of  the  expansioh 
coefficients  of  the  functions  and  Tatole  ®  contaihs 
the  coefficients. 


n 

0 

1 

2 

3 

4 

5 

6 

CM 

0 

-114.5 

-  4.8 

-33.2 

11.7 

-25.9 

3.5 

b<2> 

n 

408.5 

409.3 

-115.6 

69.6 

-15.2 

23.2 

78.3 

( 


n 

7 

8 

9 

10 

11 

b<2> 

n 

14.3 

-  0.4 

10.4 

1 .8 

-  0.6 

Table  8. 

All  values  are  to  be  multiplied  by  iC/N 
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